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ABSTRACT 



AT>£VA X 

A study of Dynamic Systems with three variable para- 
meters is made by an initial scanning of the basic geo- 
metric properties on the three dimensional space generated 
by these parameters. From these geometric properties, a 
Root-locus technique that simplifies greatly the amount of 
work in analysis and design is developed. This technique 
is extended to systems with k variable parameters. 

Finally, singular lines in the parameter plane are treated 
as a special case of the parameter space, and formulae are 
derived for a fourth order system leaving the field open 
for systems of higher order. 
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1. INTRODUCTION OF THE THIRD PARAMETER 



This study takes as a starting point Mitrovic's ap- 
proach to the root-finding problem. Although this is 
essentially a two-variable method which will yield only 
two relations out of the three needed for the three-para- 
meter case, the parametric equations obtained will help in 
understanding the nature of the basic geometric properties 
of the parameter space generated by ot, P, and y, the three 
variable parameters. The fact that it is possible to 
obtain only two relations is justified because a point on 
the s-plane is completely defined by only two conditions, 

and therefore any other possible relations will be linear- 

2 

ly related to the two originated by Mitrovic's root-find- 
ing method. 

Let us now describe briefly Mitrovic's approach, as a 
base for the present study: 

The characteristic equation for a given system can be 
written in a general form as: 



F (S) = a S n + a n S n 1 
n n-1 



+ •*»••• 4* 3...S 4* a 

1 



O 



0 



1.1 



or in a closed form as 2 



n 



F(S) 




1.2 



k=0 



If now S is expressed as: 




1.3 
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where C and are defined in figure 1.1, which is the re- 
presentation of a point in the s-plane that is defined by 

assigning values to C and u) . 

n 




Fig. 1.1 - Representation of a point 

on the s-plane 

Equation 1.3, written in another form yields: 

S = Cd(cos 0 + j sin 0) = COe^ 1.4 

—1 k 

where 0 = cos (-£) . Then S can be written: 

= oft e-^® = (cos k0 + j sin k0) , 1.5 

from which the Chebyshev functions are obvious and are: 

T, (C) = cos k0 = cos(kcos ^C) 



1.6 



£5 k = 0) k [T k (-C) + j Jl-t? U k ( -C ) ] 1.8 

= a> k [(-i) k T k (C) + j Jl-Z? (-l) k+1 u k (C) ] 1.9 

Inserting these in 1.2 and requiring that the reals and 
imaginaries go to zero independently, provides the two 
equations : 
n 

I a k“kV i)k v C) - o 

k=0 

1.1 

n 

I a k * n k (-l) k+1 VO = 0 

k=0 



But T k (o - c(u k (cj) -(Vitti) 

and upon substitution, the following equations are obtained 
n 

I<- 1)ka k"n k u k-l< ? > " 0 

k=0 

1.11 

n 

Z'- 1 > ka k"n k V f) = 0 

k=0 



where U. can be obtained from the recurrence formula : 
k 



u k+ 1 (C) = 2C u k (C) - u k . l( C) 

1.12 

where U o (C) = 0, and U^(C) =1 

Now, the coefficients are defined to be: 

a k — ^k ^k ^k 1.13 
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where a, and y are the independent variable parameters 
of the system under analysis. By substituting this ex- 
pression for a^ in equations 1.11, these equations can be 
rewritten as follows : 



+ pc 1 + yD 1 + F 1 = 0 



Q!B 2 + £c 2 + yD 2 + F 2 = 0 



1.14 



where 



n 



B 1 = h-v\ “n ku k-l (C) 



k=0 



n 



B. 



- b 



-l)\ o. n k u k (C) 



k=0 



n 



n 



C 1 Z (_1) c k W n U k-1^ ) C 2 c k ^n U k ( ^ 

k=0 



k=0 

n 



1.15 



n 



D 1 - “n ku k-l (C) E 2 - 1<- I>kd k "nX 10 



k , . k 



k=0 



k=0 



n 



■ b 



-l) k f. co k u. .(C) 

k n k-1 



k=0 



n 



- b 



-l)\ « n \<t> 



k=0 



Information concerning the design properties of the system 
must be obtained from these expressions (equations 1.14). 
This information can be: 

a). For a given set of Performance Specifications, 
what values of j3, and y, the variable parameters of the 
system, are needed in order to obtain the dominance con- 
ditions that will allow us to meet the given performance 
specifications? 
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b) . The physical meaning of the values of a, ft, and 
y, whether or not these values are actually attainable, 
and if not, what modifications must be introduced to the 
system, such that it is possible to obtain the desired 
values . 

c) . Existence of dominance conditions, if these 
exist, and the locations of the remaining roots on the 
s-plane, and their influence in the behavior of the system 

If C and 0^ are given, then the set of equations 1.14 
will contain three unknowns: Oi, f$, and y. Therefore the 
surfaces generated on the ot, ft, y coordinate system can be 
expressed only in a parametric form, generating sets of 
parameter-planes in ol and $ for a given value of y; or it 
is possible to obtain a three dimensional representation 
which is not a practical tool for engineering purposes. 
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2 . 



REAL ROOTS 



The characteristic equation for a system can be re- 
written as 

OiA(S) + j3b(s) + yc(s) + D (S) = 0 2.1 

By making S = - C^, A(S), B(S), C(S) and D(S), become real 
numbers A, B, C and D, so that the above equation becomes 

OA + 0B + yC + D = O 2.2 

which is the equation of a plane on the system of coordi- 
nates generated by a , ft, and y, with intercepts as 
follows : 

on the ol axis a = - d/A 

on the ft axis ft = - D/B 2.3 

on the y axis y = - D/C 

This is, then the first set of known surfaces that can be 
used in the three dimensional space generated by a, ft, and 
y . However, there is need for a graphical representation 
which allows the determination of the desired information 
without tedious use of computarized data, with which it is 
hard to get a clear vision of the trend of the different 
values for changes in the parameters of the system. In 
order to obtain this graphical representation, it is possi- 
ble to use descriptive geometry as a useful tool for the 
three-dimensional analysis. 

The construction rules are simple and easy to retain 
and with some practice it can be a tool that can be used 
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very often in this study. 

The first step is the representation of a plane in a 
three-dimensional space: A plane is completely defined if 

two lines belonging to such a plane are defined. It is 
possible to draw these two lines by determining the inter- 
cepts of the plane which represents S = - cr^ in the a, 
y space, with the and ay planes. (See figure 2.1a which 
is a representation of the plane S = - on the first 
octant of the coordinates axes) . Now, by rotating the y 
axis in the plane in the manner shown in figure 2.1a, 
and in amount such that the three axes lie in the plane 
generated by the axes a and P, figure 2.1b is obtained, 
which contains all the information concerning the three- 
dimensional graph. However, this time it is possible to 
work with real values and to plot all the information in 
one plane. Of course, once it is learned how to work with 
this construction method it will not be necessary to have 
the three-dimensional graphs. Let us assume we are to 
investigate whatever values of a and /3, are obtainable 
with given values of S = - cr^ and y = y^ . In the three- 
dimensional picture, we would pass a plane parallel to 
the plane aj3, through the given value of y = y^ in order 
to obtain the plane PNM. As can be seen, this is not a 
very practical method even if there exists the possibility 
of building an actual three-dimensional representation of 
the plane S = — CT^. By the construction method proposed, 
what has to be done is to project the chosen value of y 
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PLAHE 



r/ 

■v - 



-9fc 



Fig. 2.1a 



>e> 




Fig. 2.1b 

Fig. 2.1a - Representation of a plane in three dimensions 

Fig. 2.1b - Representation of a plane using projections and 
all coordinate axes in one plane 
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onto the line ab, determine the point N, and project this 
point N onto N', and finally draw a parallel line to ac 
through N'. Then the line de will be the actual represen- 
tation of the line MN of the three-dimensional graph, from 
which we can read the values of a and j3. This construction 
method can be more easily understood by the use of a simple 
example : 

Example 2.1 ; Let us assume that a system has the 
following characteristic equation 

s 4 + (3a+4/9)S 3 + (6a+3/3)S 2 + (10/9+50y)S + (5a+5y) = 0 2.4 

which can be rearranged in the form of equation 2.2 as 
follows : 

G (S) = (3S 3 +6S 2 +5) a + (4S 3 +3S 2 +10S ) & + (50S+5)y+S 4 = 0 

If S = -5 is to be investigated, then for this value of S: 

G (S) = -220a - 475/9 - 245y + 625 = 0 
and: 

625 

the OL intercept = - Z 220 = ^.84 

/-or 

the j3 intercept = - = 1.315 

and the y " = - =2.55 

Figure 2.2 is the representation of the positive axes a, £, 
and y. The negative axes can also be identified by label- 
ing the £ axis, as the -y axis, the y axis as the -/3 axis. 
Lines ab and ac are the intercepts of the plane S = -5 
with the a/9 and ay planes; dotted lines are extensions of 
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a (i orC-f) 



7 




Fig. 2.2 - S = -5 lines in the Cifiy space for Example 2.1 



18 



the lines on the negative side of the axes. Let us assume 
that we wish to investigate what values of Oi and £ are 
suitable for a given value of y, say y = -5. By setting 
y = -5 and projecting this value of y onto line ac, point 
d is obtained from which a line parallel to ab is drawn. 

This is line de, and this line contains the information 
about the values of cl and $ that will satisfy equation 2.4 
for S = -5 and y = -5. As a check, let us choose the point 
f; at this point cl = +5, P = +1.58 and y = -5, which sat- 
isfies equation 2.4 if S = -5. 

Another interesting feature of this construction me- 
thod is that due to the linearity properties, once the 
"advance" over the a axis for a given change in y is deter- 
mined, a grid can be easily drawn without the need of pro- 
jecting each value of y, by just advancing equal distances 
on the cl axis and drawing lines parallel to ab. The loci 
for y = -5.0, y = -7.5, y =+2.5 and y = 5.0 were drawn by 
this method. Similar results can be obtained by determining 
the "advance" in the 0 axis. In summary, all the necessary 
information is: the intercepts of the plane with the axes, 

the "advance" for a desired increment of one of the para- 
meters. Then the surface generated by S = -5 will be com- 
pletely described. Therefore, for a given value of S = -o^, 
the graph will look like the one in figure 2.3. When 
several values of <7 are needed, confusion can be avoided by 
identifying each pair of lines by the common point on the 
a axis (See figure 2.3). 
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Fig. 2.3 - Grid of loci for several values of real roots 
in the afiy space. 



20 



This method of construction will be useful in cases in 
which it is desired to analyze well damped systems and 
also as an easy way to construct a grid of values of real 
roots, which can be superimposed over a previously deter- 
mined set of constant C curves, in order to determine the 
possible values of real roots expected for given values of 
a, 0, y, C and 0) . 

n. 

Example 2.2 : Consider a voltage regulator. 

The compensation schemes have been worked out and this 
example will be analyzed through this study so we can ob- 
tain a better picture of the use of this construction me- 
thod. 



After manipulation of the block diagram the character- 
istic equation of the system is: 



S 4 + (2xl0 6 K 10 +1350) S 3 + (3 x10 8 K 10 +4x10 8 K i:l + 3 . 85xl0 5 ) S 2 

+ ( 10 l0 K. n +2xl0 10 K. . +3 6x10 6 ) S + 2xl0 10 K K +2.1xl0 10 = 0 
10 11 a p 

by letting K 1Q = a, = 0, K a K p = y , and rearranging the 

characteristic equation in the form of equation 2 2, it 



becomes : 

(2xl0 6 S 3 +3xl0 8 S 2 +10 10 S) a + (4x10 8 S 2 +2x10 10 S) 0 + 
2xl0 10 y+S 4 +1350S 3 +3.85xl0 5 S 2 +36xl0 6 S+2.1xl0 10 = 0 



For S = -1 the equation of the plane is : 

-9702tt - 196000 + 20,000y + 20964.38 = 0 
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where the a intercept 



20,964.38 

-9702 



2.16 



the £ intercept 



20,964.38 

-19,600 



1.06 



and the y intercept 



20,964.38 

20,000 



-1.04 



For S = -5, the equation of the plane is: 



-42750! - 9000/3 + 2000y + 2082.95 = 0 



the a intercept 




and the y intercept 



the £ intercept 




The results of this problem can be easily checked on 
figure 2.4 for the case let us take say a =0.5, £=0.5 
by drawing a line , parallel to ab through the point (a = 0.5, 
£ = 0.5), we get point c on the a axis, and then by project- 
ing this point onto the line corresponding to cr = -1.0 on 
the plane ay, we get y = -0.315 these values satisfy the 
characteristic equation of the system for S = -1.0. Another 
variation can be obtained by selecting a value of y, this 
is done for a = -5.0 in which y = 1.5 was selected, by pro- 
jecting and drawing a parallel line to de we obtain the 
line labeled (a = -5.0, y = 1.5). Any pair of values for a 
and £ on this line will satisfy the characteristic equation 
for the system. We have picked the point B, which yields 
a = 2.0 and £ = -0.385 which checks. 
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Fig. 2.4 - S = -5.0 and S = -1.0 loci for Example 2.1 
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3. COMPLEX CONJUGATE ROOTS 

Usually dominance of a certain pair of complex roots 
is desired in the system. By dominance it is understood 
that the system will behave in such a manner that its dy- 
namic performance can be approximated by a second order 
system; that is, that the influence of the remaining 
roots of the characteristic equation will be very small 
on the transient behavior as compared with the influence 
of the pair of dominant poles. This feature will allow 
the designer to obtain a certain degree of control towards 
the achievement of some performance specifications such as 
settling time, natural frequency and damping coefficient 
(overshoot or undershoot) . For this analysis equations 
1.11 can be rewritten as: 

n n 

0 !<- 1 >’Vn k U k-l< C > + * 

k=0 k=0 

n k 

+ y IwiWi'O + - 0 

k=0 k=0 

n n n 3 ‘ 1 

“ Z(-l)X“> k <0 + l> I(- 1 ) c k “n ku k (C)+ y 

k=0 k=0 k=0 

n 

+ I‘- 1 ) ke k “„ ku k <? > = 0 

k=0 

where a^ has been replaced by 

a k = ^k + ^ C k + rd k + e k* 3,2 
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By assuming given values of C and for a desired domi- 
nance condition, equations 3.1 can be written as: 



n 






n 




n 


n 


a Z b k D k 


+ 


0 


i c k°k + r 


iv k 


+ I e k D k = 0 


k=0 






k=0 




k=0 


l'.=0 


n 

r—t 






n 

r“r 




n 

< — i 


n 


“ 2-Vk 


+ 


0 


Iv* + 


y 


K E k 


+ 1 , e k E k = 0 


k=0 






k=0 




k=0 


k=0 


where 


D k 


= 


(-l) k u k 


u k. 


-1< C > 




and 


E k 


= 


(-l) k to k 
n 


u k 


(C) 





By evaluating the quantities under the summation sign these 
two equations can be further reduced to: 

0!B 1 + ^C 1 + yD 1 + E 1 = 0, and 

3.5 

0® 2 + £c 2 + yD 2 + E 2 = 0. 

These are the equations of two planes in the space OL, 0 , y. 
The locus of the values of OL, 0 and y which satisfy the 
desired values of C and C0 n is then the straight line defin- 
ing the intercept of the two planes. It is known that a 
real root can be represented as a plane in the space a 0 y, 
so now we can determine the intercept of the real-root 
plane with the straight line defined by the intercept of 
the planes described by equations 3 „ 5 . There are two 
possibilities: a) the line is not contained on the plane 

and then the intersection will be a point in the space OL 0 
y; or b) the line is contained on the plane and then we 
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have a set of values of Ot P y on the line that will N 
satisfy the desired conditions. In other words, it will 
be possible to solve a set of three equations with three 
unknowns as follows: 

From equation 3.5 

OB + ^C 1 + yD 1 + E 1 = 0 

3.6 

dB 2 + jSc 2 + yD 2 + E 2 = 0 

and from equation 2.2 

OA+j3B + yC+D = 0 
By Cramer ' s rule : 

-E 1 (C 2 C-BD 2 ) + e 2 (C 1 C-BD ) - D(C 1 D 2 -C 2 D 1 ) 

a = 5 ' 

B 1 (-E 2 C+DD 2 ) - B 2 (-E 1 C+D D 1 ) + A(-E 1 D 2 +E 2 D 1 ) 3 7 

3 = s ' 

and 

B 1 (-C 2 D+B E 2 ) - B 2 (-C 1 D+BE 1 ) + A(-C 1 E 2 +C 2 E 1 ) 

y = A ' 

where A = B 1 (C 2 C-BD 2 ) - B 2 (C 1 C-BD 1 ) + 

Let us now illustrate this use of known surfaces in the 
CL j9 y space by an examples 

Assume a third-order system with the following 
characteristic equations 

S 3 + ( 3 a +4 jS) s 2 + (4ot+3&) S + 20C + y = 0 3.8 

Assume now that we wish a pair of complex poles with 
second-order dominance at S = -5 - 8.66j and a third pole 
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at S = -50, What values of a, fi and y must we use to 
obtain the desired pole values? 

By applying S = -50.0 to equation 3 . 8 we obtain: 

7,302a + 10,000/8 - 149y - 125,000 = 0 3.9 

and by substituting S = -5 - 8.66j in equation 3 . 8 we 
obtain : 



298a + 400 j9 - y - 1,000 = 0 
and 260 a + 400j3 - 30y = 0 



3.10 



By solving the system of equations formed by 3.9 and 3.10 
we get values of 

a = 7200, j3 = -5, 385 and y = -9,400 
By substituting these values in equation 3.8, we get 
S 3 + 60S 2 + 600 + 500. = 0 

or (S + 50) (S 2 + 10S + 100) = 0 3.11 

9 

or (S + 5') (S + 5 — 8 . 66 j ) (S + 5 + 8.66j) — 0 

We see that we have three degrees of freedom and therefore 
have the opportunity of choosing three of the roots of the 
characteristic equation for systems of order higher than 
three, once we have made our choice of three of the roots. 
The remaining roots will depend on the values of a, $ and 
y (except in the very special case of situation i previously 
described, in which we can find another root by determining 
the intersection of the line with' the plane generated by 
this extra root) . Therefore, we have the freedom of choos- 
ing three of the roots of the polynomial. However, for 
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polynomials of a degree higher than three, this method does 
not tell us anything about the remaining roots of the poly- 
nomial, and once the choice has been made, the remaining 
roots will be determined. So far, as in the case of the 
parameter space, the designer has the freedom of choosing 
several of the roots of the characteristic equation, but 
is unable to determine whether or not the remaining roots 
will help in design, or even make the system stable. This 
is a tedious trial and error process, and it is very likely 
that for the chosen roots there will not be a satisfactory 
set of remaining roots. There are some useful methods like 
the Polak method and the Ross -Warren method, employing 
filter compensation, but they are essentially one-variable 
parameter methods and are not applicable in a general sense. 
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4. CONSTANT C SURFACES 



It was shown in Section 3 that the locus on the ol $ y 
space for given values of C and u> n was the line defined by 
the intersection of the two planes represented by equations 
3.5 which are: 



«B 1 + &C 1 + yD 1 + E 1 = 0 
and Q!B 2 + $C 2 + yD 2 + E 2 = 0 



3.5 



The above statement can be interpreted also as "for a given 
C, the constant loci on the ot (3 y space are straight 
lines", so that by keeping C constant for each value of 
there will be a set of two planes in the form of equations 
3.5 with the intercept a straight line in the three dimen- 
sional space Q! |S y. By projecting this straight line onto 
one of the planes generated by the coordinate axes we will 
be able to summarize all the information pertaining to this 
straight line in one plane. Figure 4.1 shows the intersec- 
tion of two planes ABC and DEJ as the line GH and its 
extension (dotted line) . The projection of GH on the plane 
is the line HF and its extension. As can be seen, equal 
increments in y will produce equal advances along the pro- 
jection. Therefore, linear interpolation can be done 
along the projection line HF, and the information concern- 
ing the line of intersection can be summarized in the a P 
plane. The line HF and its extension can easily be deter- 
mined by finding the points H and F. Setting y = 0 in 
equations 3.5, and solving for & and /9 we obtain the point 
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/ / 



I ' 

/ / 
t / 

/ 

/ 

Fig. 4.1 - Intersection of two planes in a three 
dimensional space. 
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H, applying Cramer's rule: 



(-E,C 9 + EC,) 

v - - 



and 



0H = 



^ _B 1 E 2 + B 2 E 1^ 
AH 



where 



An = Bl c 2 - 



Similarly for point F, setting $ = 0 



0!F = 



~ E 1 D 2 + E 2 D 1 
AF 



and 



0F = 0.0 



4.1 



4.3 



It is possible also, to determine a "measuring-stick" 
along the line HF by determining the value of y, given by 



( b 1 e 2 + b 2Ei ) 
y l ~ AF 



4.4 



where AF = B^D 2 = B 2 D^ 



4.5 



As can be seen, for an increment Ay = y , we must advance 
along the projection line equal distances HF. We can draw 
the projection line in the & ft plane and by using linear 
interpolation, label along this line the values of y. 
Figure 4.2 illustrates the construction method used in the 
two planes shown in the three-dimensional picture in Fi- 
gure 4.1. For this we have drawn each one of the planes 
represented in equation 3.5 by the method proposed in 
Section 2. The intersection of the two planes in the a $ 
plane for y = 0 is the point H and the intersection in the 
a /9 plane through y = y^, is the point F. Therefore, the 
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Fig. 4.2 - Two dimensional representation of the inter- 
section of the two planes shown in Fig. 4.1. 
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projection line is HF. Any other point along this line can 
be determined by: a) picking up the corresponding value of 

y^ (see y^) and projecting this value onto the a axis 
through the lines on the a y plane (ED and BA in this case) 
These projections will determine the points P and Q on the 
a axis. By drawing a line parallel to BC through P and a 
line parallel to EJ through Q we determine the point 
y = y^ = 3y 1< along the line HF, from which we can obtain 
the corresponding values of a and ft. b) Linear interpola- 
tion can be used also, by laying off the distance HF along 
the projection of the line for equal increments. Ay = y . 
The construction method was used only as a graphical help 
to understand the linear characteristics and the possibili- 
ty of linear interpolation. Everything can be done analy- 
tically. The equation of the projection line can be obtain 
ed by using the formula for a straight line, given two 
points H and F. An example follows that will illustrate 
the method. 

Example 4.1 : Assume the system shown in figure 4.3 is 

to have a dominant pair of poles with C = 0.5 and o? n = 2.0. 
The characteristic equation is 

(S + 0.5) (S + 1) (S + 5) (S + 10)+ K K S 2 + K K.S + K K. 

a L. -L 

4.6 

By rearranging terms we get: 

S 4 + 16. 5S 3 + (K K a + 73) S 2 +(K K fc + 82.5)S +(K K ± + 25)= 0 

4.7 
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Fig. 4.3 - System with three variable parameters. 
Example 4.1. 
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or 



s 4 + 16.5s 3 +(a + 73)S 2 +(j9 + 82 . 5) S + (y + 25) = 0 

where a = K K = K K. y = K K. 

3. t 



by setting C = 0.5 and = 2.0 we get the equations: 



4a - y + 135 = 0 
and 4a -2^ + 110 = 0 



4.8 



by making y = 0 we get 

and 



Q!H = -33.75 

0H = -12.0 



4.9 



For = 0 we get 

and 



ttF = -27.75 
y = 24.0 



4.10 



We can now draw the straight line on the a # plane as 
shown in figure 4.4. 

By inspection of equation 4.7 we can see that one of the 
roots will go in the right-half side of the s-plane if y 
is negative and of absolute value greater than 25; there- 
fore, for the system to be stable y must be y > -25. By 
setting y = 0 the following roots are obtained for equa- 
tion 4.6: 

2 = ~1*0 i 1.73j (desired value of roots) 

5 3 = -14.055 

5 4 = -0.446 

The corresponding values of a and are -33.75 and -12.0 
respectively. Figure 4.5 shows the location of the roots 
for several values of y. Two important features are 
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Fig. 4.4 - Straight line representation of line 

projection on the Oifi plane - Example 4.1. 
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noticeable on this graph: a) Due to the relationships 

established in equations 4, ,8, the desired pair of complex- 
conjugate roots remains in the same place on the s-plane 
for all values of y. b) The remaining roots of the poly- 
nomial travel through the s-plane in the same fashion as 
on the root-locus. This is essentially a method for deter- 
mining the loci of the roots when there is only one vari- 
able parameter in the system. This fact is justified 
because two degrees of freedom have been used in fixing 
two of the roots of the polynomial. Therefore it is pos- 
sible to obtain a root-locus for the remaining roots of 
the polynomial. The method is outlined in Section 5 and 
proves to be useful in design for systems with several 
variable parameters. Table 4.1 gives the values of the 
remaining roots and the corresponding values of oc, fi and y 
for the characteristic equation of Example 4.1. 
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Fig. 4.5 - Locus of remaining poles in characteristic 
equation of system in Example 4.1. 
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No. 


y 


S 3 


S 4 


0L 


0 


a 


0 


-14.055 


-0.446 


-33.75 


-12.0 


b 


39 


-13.29 


-1.21 


-24 


7.0 


c 


103 


-11.77 


-2.73 


- 8 


39 


d 


231 


-7 . 25+3 . 37 j 


-7 . 25-3 . 37 j 


24 


103 


e 


295 


-7 . 25+5 . 22 j 


- 7.25-5.22j 


40 


135 


f 


495 


-7 . 25+8 . 79 j 


-7 . 25-8 . 79 j 


90 


235 


g 


625 


-7. 25+10. 3j 


-7. 25-10. 3j 


122.5 


300 


h 


1007 


-7. 25+45. 3j 


-7. 25-45. 3j 


218 


491 



Table 4.1 - Corresponding values of Sy a and # for 

system of Example 4.1. 
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The straight-line property of the constant lines can 

also be used for interpolation. In order to illustrate 

3 

this let us assume that computer results yields curves 
A(C = a, y = y and B(C = a, y = y^) on figure 4.6; but 

because of some design considerations the range of admis- 

\ 

sible values for Qi and P is the region OPQR on the cn P 
plane. It is possible then to get a fairly accurate 
picture of the shape of the curves in that region by draw- 
ing straight lines between points on curves A and B having 
corresponding values of and by linear interpolation on 

these lines to obtain curves for intermediate values of y 
between y^ and y^. Curve C corresponds to y = (y^+ y^)/2 
and was obtained by dividing the straight lines into two 
equal segments. As it can be seen the loop on curve A 
tends to disappear as we approach curve B. 
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Fig. 4.6 - Use of linear interpolation in constant C 
surface using constant lines. 
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5. MODIFIED ROOT LOCUS 



There are many instances where the performance speci- 
fications for a system require a pair of complex conjugate 
roots for the system. The parameter-space techniques tell 
us what values of the parameters OL, P and y are suitable for 
the system to obtain this pair of complex conjugate roots. 
Unfortunately these techniques do not disclose any informa- 
tion about the rest of the roots of the characteristic 
equation for the given values of the parameters unless we 
solve the polynomial, and if dominance or stability condi- 
tions are not satisfied, we will have to go back and make 
another choice of values of ot P y. After several trials it 
is possible to establish the trend of the roots. However, 
this is a difficult task even if we have a computer program 
available because in general it is very hard to get a 
clear picture, and in order to get detailed results we 
would need a very complete three-dimensional picture of the 
parameter space. Because of this inconvenience we may use 
a method that gives a picture of the system at one glance 
and tells us something about the remaining roots of the 
polynomial which represents the characteristic equation. 

The name assigned to this method is the "modified root lo- 
cus", because the resulting locus follows the same tech- 
niques applied to the root locus, except that the desired 
pair of roots remains in the same place on the s-plane 
while the variable parameters Oi, fi and y are changed. 

Let us assume we have a system with characteristic 
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equation G(S, ot, y) , and a pair of complex conjugate 
roots with given values of C and C0 n is desired. It is 
possible to reduce the order of the characteristic equa- 
tion by making the division: 

g(s, <*, (L. y) 5 i 

S + 2£<j0 + CO 

n n 

The quotient of this division ♦rill be a polynomial of 
order n-2 : 

F(S, ot, 0, y) 5.2 

and the remainder will be 

F 1 (a, 0, y) s + F 2 (a, 0, y) 5.3 



The division has to be an exact one, in order for 
S + 2£o> n + to satisfy the equation G(S, Ot, ft, y) =0. 

Therefore, the remainder (Expression 5.3) has to be equal 
to zero, and for this condition to be satisfied, the values 
of F^(a, 0, y) and F 2 (0!, 0, y) have to go independently 
to zero. This can be readily seen because for any arbitrary 
value of S equation 5.3 has to go to zero. Then: 



F^ (Ot, 0, y) = 0 

and F 2 ( a » 0, 7) = 0 



5.4 



This set of equations is, in some cases, the very same set 
of equations represented in equation 3. 5, or a linear combi- 
nation of them. Therefore, they represent a line on the 
three-dimensional space Ot (3 y. We can apply the same tech- 
niques used in Section 4 to obtain a graphic or mathematic 
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representation and obtain the corresponding set of values 
of OL, 0 and y that will yield a pair of complex roots 
having values C and u.^. But what about the remaining roots 
of the polynomial? Will they yield the desired dominance 
conditions? Will the system be stable? What are the re- 
quired values of a, and y? Will these values be phy- 
sically attainable? The answer to these questions lies in 
the quotient of the division represented by Expression 5.2. 
The fact that two of the roots have been fixed means that 
we have implicitly used two degrees of freedom. Therefore, 
there is left only one degree of freedom which is the 
essence of the root-locus method. Its techniques can be 
applied using the following procedure: 

From equations 5.4, we can get any two parameters as 
a function of the third parameter. Let us assume we are to 
solve Oi and 0 as functions of y; then equations 5.4 yield: 



a = G 1 (y) 



and 



0 = G 2 (y) 



or, solving for a and y: 



and 



a = G 1 1 (0) 
y = G 2 1 (0) 



5.5 



or, solving for 0 and y: 



0 = G 1 11 (a) 



and 



y = G 2 i;L (a) . 
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Now by substituting the values of say, a and P in Expres- 
sion 5.2 (the quotient) we obtain 

F (S , y) 

or if we use a, and y: 

F(S, fi) 5.6 

or if we use $ and y: 

F (s , a) 



Now we are ready to apply the root-locus techniques to the 
polynomial F(S) . Since this polynomial has also to satis- 
fy the characteristic equation we can say that: 

F (S , y ) =0 5.7 



by rearranging terms of this polynomial we get 

r n(s) = , 

D (S ) 



5.8 



The root locus for y will be the locus of the remaining 
roots of the characteristic equation. From this root locus 
it is possible to establish stability and dominance condi- 
tions. Choose the more convenient value of y, and by sub- 
stituting it into equation 5.5 get the corresponding values 
of a and or by a simple look at the straight-line method 
established in Section 4, obtain the physically attainable 
values of a and /3. This modified root locus will also 
allow the designer to select the best choice of the remain- 
ing roots of the polynomial that will help to attain other 
performance specifications, such as bandwidth and gain. 
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The method described above for three parameters can 
be generalized for dynamic systems with k variable para- 
meters. In this case the designer has (k-1) choices of 
the roots and there will be one "free parameter" that can 
be varied at will for application of the root-locus tech- 
niques. The characteristic equation will be divided: 



G(S, 



,k-l 



a, 0, y. k) 

+ b. _S +••••• b 
k-2 o 



5.9 



where the denominator of equation 5.9 corresponds to the 
product of the chosen (k-1) roots. 

The quotient of this division will be a polynomial of 
order n - k + 1 



F(S, a, 0, y k) 



5.10 



and the remainder will be 



k)s k 2 + 



F k-2 (a ' ^ ,y k)S + F k-l (fl ^' y,,,) 

5.11 



Again, if the division is going to be satisfied with the 
remainder equal to zero a set of k-1 simultaneous equa- 
tions are obtained. 



F-^a, 0, y ... k) = 0 
F 2 (a, 0, y ... k) = 0 



F k-2 (Qi ' 7 “• k) = 0 

F k_i (a, 0, y ... k) = 0 



5.12 
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From equations 5.12, it is possible to solve for (k-1) para- 

til. 

meters as functions of a k parameter and by substitution 
in Expression 5.10 we get: 

F (S , free parameter) =0 5.13 

The application of root locus techniques to equation 5.13 
will produce the locus of the remaining roots of the cha- 
racteristic equation. These techniques are useful in va- 
rious cases, depending on the applications: 

a) . for two variable parameters it is possible to set 
one of the roots of the polynomial at a large negative value 
and obtain a root locus for the remaining roots of the po- 
lynomial. 

b) . for three parameters the choice of the dominant 
pair of complex conjugate roots can be made, allowing de- 
termination of the remaining roots of the polynomial. It 
is possible to give two of the real roots large negative 
values and determine the remaining roots from the root locus, 
choose the ones that help us to get as close as possible to 
the desired performance specifications. Similar considera- 
tions can be made for k variable parameters. 

The calculations involved are simple and provide an 
overall picture of the system and its possibilities. In 
many cases the quotient is a function of only one of the 
parameters; therefore, the equation for the root locus is 
available without the need of substitutions (see Problem 
5.2). With a computer program for root locus the job is 
less laborious. 
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Several examples are available in order to clarify the 
procedure and the techniques. Each one of the cases out- 
lined above is studied. 

Example 5.1 ; 




System with three variable parameters - Example 5.1. 

Fig. 5.1 

The system on Figure 5.1 is to be compensated using tacho- 
meter feedback and accelerometer feedback. We wish to in- 
sure dominance of a pair of complex conjugate roots with 

C = 0.5 and <j0 =2.0. We wish also to know what are the 

n 

permissible values of the remaining roots of the character- 
istic equations and, if it is possible, to obtain dominance 
conditions the required values of p, k and k. . The trans- 
fer function for the system on Figure 5.1 is: 

FC3) = 100 

S 5 +(09-50)#+(50a+830)S 3 + (830a+100£+6250)S 2 + (€250a+100yf 1 500)&KL5300aH-2 5 ) 

5.9 
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where p = a, K = 0, and K. = y 

d t 

. . 2 

By dividing by S + 2S + 4 we get the quotient: 

S 3 + (a+48)S 2 + (480! + 7 30) S + (7300! + 100/3 + 4598) 

5.10 

and the remainder 

(4598G!+100y - 200 £+2884) S + 120800! - 4000 - 18292 

5.11 

In order for the remainder to be zero for arbitrary values 
of S, it is necessary that 

A 

45980! + 100y - 2000 + 2884 = 0 

5.12 

and 1208tt - 4000 - 18292 = 0 

So far, two degrees of freedom have been used; therefore 
we are left with only one of the parameters variable (free 
parameter) . The other two are dependent functions of this 
variable parameter (equations 5.12). 

It is easiest to have a as the variable parameter, 
because most of the coefficients in the quotient, Expres- 
sion 5.11, are functions of a, and also because 0 is 
easily expressed as a function of a, in the second of 
equations 5.12. From these considerations: 

o 3020a - 4573 r „ 

0 = jQo 5.13 

By substituting the value of 0 in Expression 5.10 and 
making it equal to zero, we get: 

S 3 +(a+48)S 2 + (48+7 30) S + ( 3750 a + 25) = 0 5.14 
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Now, by applying the root locus technique, 

a(S 2 +48S+3750) =_ i 

S 3 +48S 2 +730S+25 

a(S + 24 - 56 . 338j ) (S + 24 + 56.338-i) 

° r (S+23 . 847-jl2.377) (S+23 . 847+j 12 . 377) (S+0.0343) 1 

The root locus for the variable parameter Oi is shown in 
figure 5.2. For any point on the given root locus it is 
possible to determine the required value of Oi and the cor- 
responding values of & and y are determined from the rela- 
tions expressed in equation 5.12. The graphical method 
may also be used. The root locus in figure 5.2 shows all 
possible values of the remaining three roots of the poly- 
nomial. For low values of Oi, the real root will be very 
close to the pair of complex poles that are desired as do- 
minant. For high values of Oi, the real root will go far 
in the left of the s-plane and the other pair of complex 
roots will be very close to the generating zero of the 
root locus. As can be seen, the system allows in this case 
an improvement in the bandwidth by choosing the value of ol 
that will make the real root lie in the far left of the 
s-plane and selection of the complex roots to insure a 
dominance condition (real part of the complex roots about 
ten times the real part of the dominant roots) . Other 
considerations can be made such as what values of Oi will 
yield permissible values of K 2 and K fc , without saturating 
the amplifiers that will provide these gains, etc. 
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Fig. 5.2 - Modified root locus for remaining roots of 

characteristic equation for system in Fig. 5.1. 
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Example 5.2 ; 

A system has the following characteristic equation: 

G (S) = S 5 + 16S 4 + 40S 3 + (a+10) S 2 + ( $+50) S + (7+25) 
It is desired to have a pair of complex conjugate roots 
with C =0.5 and to = 2.0. By making the division: 



s 5 + 16S 4 + 40S 3 + (a+10) s 2 + ( $+50) S +7+25 

S3 + 2S + 4 

we obtain the quotient: 


5.15 


S 3 + 14S 2 + 8S + (a-62) 
and the remainder 


5.16 


(-2a+0+142) s + (-4^+7+273) = 0 
By making the remainder equal to zero we get: 


5.17 


- 2 a + 0 + 142 = 0 
and 

- 4 a + 7 + 273 = 0 


5.18 



As is easy to see in this case there is no need for substi- 
tution, since the quotient is a function only of a. However, 
the values of j? and y will depend on the value chosen for a, 
from the relationships expressed in equations 5.18. 
Rearranging Expression 5.16 for root locus application: 



S 3 + 14S 2 + 8 S + (a-62) = 0 


5.19 


or 

* = - i 

S 3 + 14S 2 + 8 S - 62 


5.20 


or 

a - = 1 
(S - 1.74646) (S + 2.7266) (S + 13.01977) 


5.21 
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The root locus for equation 5.21 is shown in figure 5.3. 

As can be seen it is not possible to obtain dominance con- 
ditions. The system will become unstable for very large 
or very small values of ol. in this case, with very few 
steps it was possible to obtain a clear picture of the 
system. It is obvious that the root locus must have the 
same shape whatever parameter is chosen. We can check 
this fact by choosing say, y as our variable parameter. 
From equations 5.18 the expression for ol is: 



<* = 2L±J2I 5 . 

By substitution of the expression for ol in the quotient 
(Equation 5.19) we get 

4S 3 + 56S 2 + 32S + y + 273 - 248 = 0 
or 

Z = _ 1 

4(S 3 + 14S 2 + 8S + 6.25 



If K = y/4, we get: 

K 

(S 3 + 14S 2 + 8S + 6.25) 



1 



or 



(S + 13.439) (S + 0.28 + j0.622) 



5.23 



The root locus for equation 5.23 is shown in figure 5.4 and 
it is exactly the same as the one on figure 5.3. The part 
which is not present is that corresponding to negative 
values of y. Had we drawn the root locus for negative 
values of y this part would appear. The points where the 
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Fig. 5.3 - Modified root locus for system in Example 
Q! as free parameter. 




Fig. 5.4 - Modified root locus for system in Example 5.3 - y as free parameter. 




UJ 

< 

\ 

<T) 



-r 

1 l 



4 GO 
I 



i 

4 



04 







i 



55 



-/o 



generating poles for the root locus of a appear are the 
points for which the value of y make a = 0. Similarly 
the points where the generating poles for the root locus 
of y appear correspond to the value of ol that makes y = 0. 
We can make a quick check. Let us determine the value of 
at at point A on figure 5.3. This value yields a = 68.25, 
and by substituting in equation 5.22 we find that this is 
exactly the value of a that makes y = 0. 

Example 5.3 ; 

Let us consider again the example 4.1 in which the charac- 
teristic equation is 

S 4 + 16. 5S 3 + (a+73)S 2 + (j9 +82 . 5) S + (y+25) = 0 

by using the same desired values of C = 0.5 and = 2.0 

2 

and dividing by S +2S+4, the quotient is 
S 2 + 14. 5S + (a+40) 
and the remainder is 

(-20!+ 0-55. 5) S + (-4a+y-135) 

Setting the remainder equal to zero, we have 
-40! + y -135 = 0 

, 5.24 

and 

-2oi + 0 - 55.5 = 0 

By comparing equations 5.24 with equations 4.8, in example 

4.1, it is easy to see that these two sets of equations are 

exactly the same; therefore, equations 5.24 represent the 

straight line that defines a constant u) on the constant C 

n 
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surfaces. By making a root locus set up with a as a free 
parameter we get: 



(S 2 +14 . 5+40) 



(S+10.995) (S+3.705) 



5.25 



The root locus is shown as the slashed lines in figure 5.5. 
Poles marked D are the generating poles using a as a free 
parameter. Had y been chosen as the free parameter the 
generating poles would be the ones marked as B on figure 
5.5, and the root locus using y as the free parameter would 
be the part shown as heavy lines, plus the slashed lines. 
For roots at points between B and D the value of & would be 
negative. This fact can be checked in figure 4.4 where the 
corresponding points B and D are marked. The point of im- 
portance in this case is that since root locus techniques 
consider only positive values of the variable parameter, 
for each parameter chosen as free parameters, we will get 
only part of the overall root locus. Therefore, if the 
results with a positive value of the free parameter are not 
satisfactory, it is good practice to make a quick check for 
negative values of the free parameter. This can be done 
easily by starting at the generating poles and going in the 
opposite direction. For this specific example, we can see 
that for negative values of the root locus will go from 
D towards B(see Figure 5.3), and eventually one of the roots 
will cross the j 0) axis and the system will become unstable. 
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Fig. 5.5 - Modified root locus for Example 5.3. 
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Example 5.4- Example with parameters product . 

There are some cases where the characteristic equation con- 
tains products of parameters. This technique is applicable 
only if the quotient of the division is a function of only 
one of the parameters. Let us assume a system with charac- 
teristic equation: 

s 4 + (a+15)S 3 + (20a+200)S 2 + (y+600) S + a/9 + 6000 

5.26 

and the system is required to have a pair of complex con- 
jugate roots with C =0.5 and = 10 . From the division 
we get quotient 

S 2 + (a+5)S + (10a+50) 5.27 

and the remainder 

(-200Q!+y-400) S + (a/3-loOOa+lOOO) =0 5.28 

The root locus technique can be applied to the Expression 

5.27, thus, 

tt(S+10) = _ i 

(S 2 +5S+50) 

or 

a(s+io) = _ i 

(S+2 . 5-6 . 644 j ) (S+2 . 5+6 . 644 j ) 
and the root locus can be seen in figure 5.6 
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Fig. 5.6 - Modified root locus for Example 5.5 - Oi as free 
parameter. 
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The choice for dominance can be made as a pair of poles at 
S = -20. From the root locus we get 

a = 35 

and from Expression 5.28: 



and 




1000a - 1000 
a 

34000 = 9?1 



y = 200 a + 400 = 7400 



we could also obtain the values by multiplication. 

G (S) + (S+20) 2 (S 2 +10S+100) = 0 

5 . 2 

G ( S) = S 4 +50S 3 + 900S 2 + 8000S+40 . 000 = 0 

By comparison of the two expressions for G(S) (Eqs. 5.26 
and 5.29) we get 



a + 15 = 50 



or 




• 


O' = 35 




20a + 200 


-- 900 




and 


y + 600 = 


8000 : 


y = 7,400 


and 


Q!j9 + 6000 


= 40.000 


: 0 = 34 °°° = 971 



which checks with the results above. 

This is a special case where the quotient (Eq. 5.27) 

is a function of only one of the parameters. In other 

cases the root locus technique is not applicable because 

2 3 

the free parameter appears in the quotient as K, K , K , 
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etc., depending on the complexity of the parameter products 
in the characteristic equation. However, there exists the 
possibility that the designer has enough freedom to modify 
the system configuration and make the selection of the va- 
riable parameters in such a way that no products appear, 
or in case of products, that the quotient be a function of 
only one parameter. All this can be achieved by changing 
some feedback paths or series compensation schemes. 
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C. SINGULAR LINES AND SELF-ADAPTIVE SYSTEMS 



In the analysis of systems with two variable parameters 
by application of the transformations proposed by Mitrovic 
(coefficient plane) and Siljak (parameter plane) , the pre- 
sence of systems with what is known as singular lines was 
detected. A summary of the procedure that led to the 
detection of this type of system follows. This leads to a 
better understanding of what is meant by singular lines. 

From Section 1, the characteristic equation for a given 
system can be written as: 

F(S) = a S n + a , S n-1 + + a,S + a =0 

n n-1 1 o 

6.1 

Repeating the manipulations of the above equation outlined 
in Section 1, the following expressions are obtained 
n 

U k-l <C) * 0 

k=0 6 . 2 
n 

I(- 1 > ka k"n k V C) “ 0 

k=0 

Now by defining the coefficients as 



a k = ^k + ^ C k + \ 6,3 

and the substitution of a^ in Equation 6.2, this set of 
equations can be rewritten as the system of linear equa- 
tions 

hB 1 + 0c x + Dj^ = 0 

6.4 

aB 2 + 0c 2 + d 2 = 0 
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where C^, C 2 , and D 2 are defined as in Equa- 

tions 1.15. When values of C and 0)^ are assigned, and the 
values of B^, B 2 » C^, C 2 , and D 2 are calculated, each 
of the equations 6.4 will represent a straight line on the 
Oifi plane and in a geometric sense, the intersection of 
these two straight lines will yield the values of a and j9 
that satisfy the system of equations 6.4, and insure that 
the characteristic equation will have a pair of complex con- 
jugate roots with the assigned values of C and But the 

intersection of two lines in a plane is not always a point. 
Two other cases are possible: a) the two lines are parallel 

in the a/3 plane and in this case no pair of a, values will 
satisfy the system of equations 6.4 (see Fig. 6.1b), or 
b) the two lines are the same in which case they intercept 
on an infinite number of points and any point along this 
line will yield a pair of complex conjugate roots with va- 
lues C and cc . The three different cases are shown in Fi- 
n 

gure 6.1. Figure 6.1c shows the case of coincidence of the 

two lines; this is known as a singular line case and its 

4 

name is derived from the matrix algebra manipulation of 
the system of equations 6.4. The systems with singular line 
characteristics find application in self-adaptive systems 
where changes in the dynamic behavior depends on the ability 
of the system to cope with variations of the parameters. 

The identification of the values of the parameters and its 
changes in these self-adaptive systems is done by observa- 
tion of the output of the system. Identification is a pro- 
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4A 




Fig. 6.1a 





Different cases of mapping of a point on the s-plane 
in the aft plane. 
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blem beyond the scope of this work, but literature on the 
subject is available (see reference 5 and 6) . Once the 
change on the parameters is identified the self-adaptive 
systems use sensitivity models which are fairly good for 
very small changes in the parameters. This is because of 
the approximations made in the selection of the sensitivity 
models. The idea of using systems with singular lines in a 
self-sdaptive system came up because it yields a precise 
mathematical expression (equation of a straight line) that 
relates the changing parameters. This expression can be 
simulated by analog or digital means, so the problem of 
steadyness of the dynamic behavior of the system is solved 
in part. A sensitivity study on the remaining roots of 
the polynomial will tell if the desired conditions are 
going to be jeopardized by the given changes in the para- 
meters and what tolerances can be allowed in the change of 
the parameters. The attractiveness of the use of systems 
with singular lines in self-adaptive systems lies in the 
fact that it is only necessary that one of the parameters 
be adjusted for changes detected in the other parameter 
and vice-versa. As was said before, the expression regula- 
ting these adjustments is the equation of a straight line, 
and if sensitivity studies are feasible, the tolerance in 
the changes of the parameters can be determined. 

So far, the detection of systems with singular line 
characteristics is a difficult job, and even in the case 
where singular lines are obtained the values of C and 0^ 
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which generate these singular lines might not be satisfac- 
tory for the desired performance specifications. The sub- 
ject is relatively new and research is going on on the 
matter. Therefore, the possibility of an easy method of 
detection and set up of singular lines might appear with 

time. There are computer programs available for detection 

4 

of singular lines . 

The results obtained in Section 4 and 5 of this work 
led to the thought that singular lines on systems with two 
variable parameters can be considered as special cases of 
a system with three variable parameters in which one of the 
parameters is independent of the other two and has a fixed 
value. This value is obtained from one of the planes re- 
presented by equations 3.5. It is clear that the above men- 
tioned plane is parallel to the one generated by the axes 
represented by the other two parameters. 

Figures 6.2 and 6.3 will help to clarify the concept 
above stated. Figure 6.2 shows the general case where £ 
and have been given and the point corresponding to the 
given value of S maps in the three-dimensional space as 
the line AB which is the locus of the values of a, P and y 
that will produce the desired pair of roots with parameters 
for S : £ and on the s-plane. Line AC is the projection 
of line AB onto the a/3 plane. Figure 6.3 shows the special 
case where plane A is parallel to the a$ plane through y = y^. 
The intercept of planes A and B is the line AB, that esta- 
blishes the relation between a and /3 as a straight line on 
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Fig. 6.2 - General case of mapping of a point on the 
s-plane into a straight line in the three 
dimensional space aft y. 
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6.3 



Special case of mapping of a point on the 
s-plane with generation of singular line. 
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the oc ft plane that passes through y = y^. (Plane A in this 
case) . So much for the geometric interpretation; the ma- 
thematical interpretation in terms of the division method 
outlined in Section 5 is that from one of the two equations 
obtained from the remainder, F^(a, P, y) = 0 or F 2 (a, 0,y) = 0 
a specific value of one of the parameters is obtained and 
by substitution of this value in the remaining equation, we 
will get the equation of the straight line relating to the 
other two parameters. An example follows. 

Example 6.1 : 

Let us assume we have a system with characteristic equation 
G ( S) = S 4 +(lOa+lO)S 3 +(4Oa+50+lOy+3O)S 2 +(8Oa+lOj3+252) 

+ (70a+10j3+25y+340) 6.5 

and also assume that somehow we have found that a pair of 

complex conjugate roots with C = 0.5 and = 2.0 will 

yield equation F^(a,j9,y) = 0 from the set of equations 5.4, 

defining a plane parallel to the a P plane. Division of 

2 

equation 6.5 by S + 2S + 4 yields as quotient: 



S 2 + (10a + 8) S + (20a + 5)9 + lOy + 10) 6.6 

and the remainder is 

(-20y + 200) S + (-10a - io/3 - 15y + 300) 6.7 

By letting the remainder equal to zero we get 

y - 10 = 0 6.8a 

10 a + 10/3 + 15y - 300 = 0 6.8b 



The graphical representation of equations 6.8 is shown in 
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figure 6.4. Now, by continuing, from equation 6.8a: 

y = 10 6.9a 

and from equation 6.8b, substituting this value of y: 

a + 0 - 15 = 0 6.9b 

Now substituting back in the characteristic equation the 
value of y we get 

(^(S) = S 4 +(lOa+lO)S 3 +(4Oa+50+13O) 3^(80^+100+ 252)S + (70tt+10£f 590) 

e.io 

2 

and dividing by S + 2S + 4 we get the quotient 

S 2 + (10a+8)S + (20a + 50 +110) 6.11 

Expression 6.11 is the same as 6.6, had the value of y been 
substituted back in this expression directly. The remainder 
for the above division is: 

-100! - 100 +150 

and by making the remainder equal to zero : 

a + 0 - 15 = 0 6.12 

This is a singular line case in the 0! 0 plane, because by 
assigning to y the value of 10 there are only two degrees 
of freedom left: 0! and 0 in equation 6.10. According to 
what was stated in Section 5 and taking G^(S) of equation 
6.10 as the characteristic equation of the system with two 
variable parameters and using the two degrees of freedom in 
setting the value of a pair of complex conjugate roots 
S = -1.0 - 1.732j, the remainder of the division should be 
a set of two linear simultaneous equations with a unique 
solution for a and 0. 
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Fig. 6.4 - Singular line case for system represented by 
Equation 6.5 - Example 6.1. 
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Equation 6.12 is the only requirement for the remainder to 
be zero and therefore an infinite number of sets of values 
of a and j9 can be assigned. The root locus for the remain- 
ing roots of the polynomial is shown in figure 6.5, with a 
as the free parameter. This was obtained by substitution of 

P = 15 - Oi (Equation 6.9b) 

into equation 6.11; then: 

S 2 + (10a + 8) S + (15a + 185) = 0 6.13 

or 

10a(S + 1.5) 1 

(S 2 + 8S+ 185) 

If 10a = k then 



K ( S+l .5) _ . 

(S+4-jl3) (S+4+jl3) 

By inspection of the root locus, and with the criterion of 

existence of dominance being that the real part of all of 

the remaining roots of the polynomial must be greater or 

2 

equal to 10 times the real part of S + 2S + 4, it can be 
seen that the root locus enters the desired region at 
point A and a section of it leaves the region at point B. 
Therefore, for values of a between these two regions we 
have acceptable dominance conditions. For point A the 
value of a is 2.41. By setting a: 



a 



OB + OA 
2 



3.61 

2 



1.80 
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Fig. 6.5 - Modified root locus a as free parameter. 
Example 6.1 - Singular line case. 
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we can have tolerances for ft of + 33% from the chosen 
value. Then for 0 we get: 

0 = 13.19 

with tolerance for 0 of - 4.65% from the assigned value. 

It has been proven that the singular line in the a 0 
plane can be taken as a special case of the three-parame- 
ter case. Unfortunately, the detection of planes parallel 
to either of the planes containing the coordinate axes is 
a 30 b as difficult as the detection of singular line in the 
Oi 0 case. From the analytic point of view it will be pos- 
sible to modify the coefficients of the different powers of 
S in such a way that a plane originally at an angle with 
the planes containing the coordinate axes becomes parallel 
to one of them. But from the practical point of view it is 
difficult to get this modification into the system, because 
these coefficients come from manipulation on the block dia- 
gram and usually affect several terms of the characteristic 
equation. Another limitation may be that the desired values 
are not physically attainable, due to the intrinsic compo- 
sition of the elements of the system. From preceeding sec- 
tions it has been known that a point on the s -plane maps 
in a straight line in the a 0 y space and this straight line 
has a projection on the ot 0 plane. Therefore, it might be 
easier from the designer's point of view, to "simulate" a 
singular line case by using three variable parameters, and 
instead of one set-up for adjustment of parameters, use two; 
that is; for each change in say, y, it would be necessary 
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to modify accordingly a and j3. The expressions that will 
determine the modification procedure are precise mathemati- 
cal expressions represented by equations 5.4 and they can 
be simulated by analog or digital means available. How- 
ever, there are considerations such as weight, space and 
economy that will put a burden on a design of this type, 
but in the event that they can be overlooked, we will gain 
the ability of choosing the pair of dominant poles that 
will help in the achievement of the performance specifica- 
tions and obtaining solid tolerance limits on the parame- 
ters. Engineering judgement must be used in this case to 
weight all these considerations. 

When the design considerations stated above become a 
real problem and a singular line case is desirable, it is 
worth while to study the coefficients of the polynomial 
which represents the characteristic equation, and by alge- 
braic manipulation, obtain the singular line case. The 
method outlined in the following section for a fourth order 
polynomial can be used in higher-order polynomials. The 
results obtained are quite satisfactory and a clear picture 
of what has to be done to obtain singular lines is presented. 
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6.1 GENERATION OF SINGULAR LINES IN A FOURTH-ORDER 
POLYNOMIAL WITH THREE VARIABLE PARAMETERS. 

A Fourth-Order characteristic equation of a system 
with three variable parameters can be written in general 
form as : 

S + (b^tt+c^ $+e 3 y+f 3 ) S 2 + (b 2 /3+c 2 $+e 2 y+f 2 ) 

6.14 

(b, a+c n 0+e. y+f ) S + (b oc+c 0+e y+f ) = 0 
1 1 1' ' K o o o' o' 

By making the substitution of 

^ = e k y + f k 6 - 15 

equation 6.14 can be written as: 

S^ + (b^tt+c ^ ^+d^ ) S^ + (b^tt+c 2 £+d 2 ) ® + (b^a+c^$+d^) S + 

(b a+c 0+d ) = 0 6,16 

It is clear that the case of singular lines will depend on 
the configuration of the polynomial, that is, on the coeffi- 
cients of the different powers of S. Therefore, it is pos- 
sible to determine what coefficients must be affected in 
order to obtain singular lines. If it is assumed that the 

polynomial has singular lines for values of x = 2£o^ and 

2 2 
y = , division of the polynomial by: S + x S + y yields 

the follwoing results: 

Quotient : 

S 2 + (b 3 a+c 3 j3+d 3 ~x) S + (b 2 -xb 3 )a +(c 2 ~xc 3 )£ + (d 2 ~y-xd 3 +x 2 ) 

6.17 
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The remainder is : 



F^at'fts + F 2 (a,fi) 



6.18 



where : 




6.19 



and 



F 2 (a, P) = (b o -yb 2 +xyb 3 ) a + (c Q -yc 2 +xyc 3 ) 0 + 




6.20 



Singular lines can be generated for each of the following 
cases : 

Case I : F-^(a,/3) = 0 for arbitrary values of 0i and /3. Then: 



which is a system with three equations and two unknowns (x 
and y) . From inspection of the system the following conclu- 
sions are drawn: 

a) . It is possible to solve for x and y by utilizing 
only equations 6.21 and 6.22. 

b) . Equation 6.23 is a function only of d^, d 2 and 
d 3 » which in turn are functions of y (equation 6.15) . 
Therefore, once the values of x and y are found, y can be 
set at the value that will satisfy equation 6.23 for the 



2 

b^ - yb 3 - xb 2 + x b 3 = 0 



6.21 




6.22 



d^ - yd 3 + 2xy - xd 2 + x^d 3 - x^ = 0 



6.23 
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values obtained for x and y. It is worth noting that the 
modifications in equation 6.23 do not affect the values of 
x and y obtained before. Then, from equations 6.21 and 
6 . 22 : 



2£w = 



n 



x = 



to 

n 



= y = 



C l b 3 - b l C 3 
C 2 b 3 “ b 2 C 3 

b ! " xb 2 + x2b 3 



6.24 



and the third equation to be satisfied once x and y have 
been found is : 

d^ - yd^ + 2xy - xd 3 + x^d^ - x^ = 0 6.25 



By substituting d^ = e^y + f^ and using the values obtained 

for x and y we can solve for y. 

If the values of C and to are not satisfactory to the 

designer, it is possible to know exactly what coefficients 

can be modified in order to get the desired values. In this 

2 

case if for example y = to^ is satisfactory but we wish to 
change x = 2^0^ , then C^, C 2 , b^ and C 3 can be modified at 
will without disturbing the value of y. From the above 
coefficients the designer might be able to choose the most 
convenient methods for the acutal modification, such as 
feedback paths, series compensation devices, etc. The 
above considerations are applicable to cases II and III 
which follow. 
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Case II. 



F 0 [Oi,fi) = 0 for arbitrary values of & and j9. Then 



b o ~ yb 2 + Xyb 3 = 0 



c o " Y c 2 + x Y c 3 = 0 

2 2 

d Q - yd 2 + y + xyd 3 - yx =0 



From the first two equations we get: 



6.26 



b c 0 - b 0 c 

20 = x = r-2-2 

n b c_ - b_c 

o 3 3 o 

~ b 

, . 2 o 

tO = y = — — — 

n J b„ - xb„ 



6.27 



and 



2 2 
d o " yd 2 + Y + *yd 3 " Y x 



= 0 



6.28 



Again, by substitution of d^ = e^y + f^, and using the 
values obtained for x and y, it is possible to solve for y. 
Case III . F 1 (a, j3) = F 2 (a, 0) for arbitrary values of a 
and 0. Then: 



b^ - yb 3 - xb 2 + x b 3 = b Q - yb 2 + xyb 3 

, 2 

C 1 -yC 3 -xc 2 +XC 3 =C o -yC 2 + xyC 3 



6.29 



2 3 2 2 

d^ - yd 3 + 2xy - xd 2 4 x d 3 - x = d Q -yd 2 4 y 4 xyd 3 ~yx 



From the first two equations we get 



_ (b 2 - b 3 ) (c 1 - c o ) 4- (b x - b o ) (c 3 - c 2 ) 

n x b 3 (c^ - c 0 - c^) 4 c^(b 0 - b n 4- b^) 



CO = y = 

n 1 



b. - xb„ 4 x b_ - b 
1 2 3 o 

b 3 + xb 3 - b 2 



'3 ' 2 



6.30 
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and 



2 3 2 2 

- yd^ + 2xy - xd 2 + x d^ - x = d Q - yd 2 + y + xyd^- yx 

6.31 

That with the given values of x and y will yield the value 
of y. 

The important feature of the results stated above is 
that the designer knows exactly what coefficients of the 
characteristic equation affect directly the values of C 
and and hence, what parts of the actual system must be 
modified in order to obtain the desired values. 

The solutions proposed above are for the case where 
the "free parameter" is y ; that is: one of the planes 

considered in Section 6 is parallel to the Oil 3 plane through 
the calculated value of y. The quotient (equation 6.17) 
will be then a function of Oi and ft. From the remainder, Oi 
can be expressed as a function of P, or ft as function of Oi, 
and root locus techniques applied in order to determine the 
location of the remaining roots of the polynomial. A sensi- 
tivity study can be made on the root locus. For applica- 
tions in self-adaptive systems there will be need of only 
one set-up for modifications when changes in either of the 
parameters occur. The following example will illustrate 
the technique. 

Example 6.2 ; 

A system has the following characteristic equation: 

S 4 + (10a+10)S 3 + (400^+5^+130) S 2 + (80O'+10i3+10y+5) S + 



6.32 
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(70O!+10£+y+200) 



Question; What values of C and c o are attainable in sin- 

n 

gular line cases and what are the required values of y? 

For Case I, substituting the values of the coefficients 
in equations 6.24; 

x = 2£co =2 or C = 0.5 
n 

2 

y=C0=4orC0=2 6.33 

J n n 

y = 24.7 

Substituting the value of y in equation 6.32 and dividing 
2 

by S + 2S + 4 we get: 

Quotient : 

S 2 + (10a+8)S + (200H-5/9+110) = 0 6.34 

Remainder; 

- 10a - 10/3 - 215.3 = 0 6.35 

Now root locus techniques can be applied to equation 6.34 
after substituting ft as a function of a from equation 6.35. 
The expressions for the quotient and the remainder can also 
be obtained from equations 6.17 and 6.20 simply by substi- 
tuting the values of the coefficients. 

For Case II, substituting the values of the coefficients 
in equations 6.27; 



x = 2CiC 
n 


= 0.5 


or C = 0.1765 




II 

to 


= 0.2 


or to = 1.4142 
n 


6.36 



y = 46.5 



Substituting the value of the y in equation 6.32 and divid- 
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2 

ing by S + 0„5S + 2 we get: 

Quotient : 

S 2 + (100+9. 5) S + (350+50+123. 25) = 0 
Remainder : 

(42.50 + 7 . 5 £ + 389. 375) S = 0 

For Case III, substituting the values of the 
in equations 6.30: 

2£to = x = 1 or C = 0.5 
n 

2 

to = y = 1 or to - 1 
n 2 n 

y = 204/9 = 22.667 

Division yields 
Quotient : 

S 2 + (100+9) S + (300+5 $+120 ) = 0 
Remainder : 

(400+50+102.667) S + (400+50+102.667) = 



6.37 



coefficients 



6.38 



6.39 



6.40 
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7. CONCLUSIONS AND RECOMMENDATIONS FOR FURTHER WORK 

In this work it is shown that the increase in degrees 
of freedom in a system really helps the designer in the 
task of obtaining from it the desired characteristics. 
Mitrovic's approach has proved to be useful in the parame- 
ter plane or when products of parameters are present in 
the characteristic equation of a system. However, when the 
number of variable parameters is three or more it is diffi- 
cult to obtain graphical information. The tabulated data 
obtainable from a computer comes only in quantized form and 
it is difficult for the designer to use it in seeking his 
objective, which is to make the best choice of the variable 
parameters. In the field of singular lines a method of ge- 
neration of singular lines for a fourth-order polynomial 
was established and it is expected that formulae can be 
obtained for systems with characteristic equation of higher 

order. With this method the values of C and 0) at which 

n 

singular lines are expected and the value at which a third 
parameter must be set are obtained. In addition, the mo- 
difications in the system for which the values of C and 
are changed are indicated. 
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